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QCD point-to-point correlation functions at distances .2-1 fm are very different for different channels, and 
they tell us a lot about inter-quark interactions. Recent studies based on experimental data, 'instanton liquid' 
approach and lattice measurements are reviewed. Agreement between all of them show that instanton-induced 
forces dominate the light quark physics, and new findings by Negele et al, that hadrons survive 'cooling', make 
this statement obvious. We also argue that chiral symmetery restoration is due to breaking of the 'instanton 
liquid' into polarized 'instanton molecules'. 



1. INTRODUCTION 

Although space is limited, let me start with a 
very general comments. There is no need to re- 
peat at this conference how complicated is the 
QCD vacuum, and how basic it is for the whole 
domain of particle and nuclear physics. How- 
ever, it seems worth repeating that lattice simula- 
tions should not only try to reproduce the Particle 
Data Table, with the best accuracy possible, they 
should lead to understanding of its structure. 

How one can do it? Of course, there are mul- 
tiple examples in other areas. Just the preceed- 
ing talk gives an example of similar attemts: cos- 
mologists try hard to find 'large scale structures' 
in positions of galaxies and inhomogeneous back- 
ground radiation. Much more familiar example 
is meteorology: very complicated charts of wind 
velocity, pressure and temperature measured in 
many points can be translated to easily unde- 
standable forecasts, just by pointing out few main 
cyclons, their magnitude and direction of motion. 

Similarly, one may look for 'structures' among 
variuos fluctuations of the gauge and quark fields, 
trying to single out the most important ones. 
I am going to argue below, that it is the 4- 
dimensional relatives of cyclons, the instantons, 
which one should look for first. 

My second remark deals with the tools used, 
the point-to-point correlation functions. We need 
them because for light quark channels (unlike 
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the heavy quarkonia) the knowledge of hadronic 
masses does not provide simple and unique pic- 
ture of interquark interactions. (Note here anal- 
ogy with nuclear physics: compare 30's, when 
only binding of nuclei were known, to the days af- 
ter precise NN scattering experiments were done.) 

Lattice simulations have an important role 
here: new mass measurements should be suppli- 
mented by new wide scale studies of details un- 
available in 'real' experiments, such as correla- 
tors, wave functions and structure functions for 
channels with various quantum numbers. 

Of particular interest is a transition region, 
X ~ 0.2 — 0.5 fm where free propagation of quarks 
and gluons (at small distances) turns into com- 
plicated non-perturbative behaviour. Discussion 
of relevent physical questions and available phe- 
nomenological information on this region can be 
found in my recent review j^]. 

Tunneling phenomena in gauge theories, dis- 
covered in , fascinating semiclassical theory, ex- 
planation of chiral anomalies |3| , first applications 
to QCD problems Q| etc. have attracted a lot of 
attention in late 70 's. However, as no explana- 
tion for 'diluteness' and validity of semiclassical 
approximation were suggested from the first prin- 
ciples, optimism has soon died out and most peo- 
ple left the field. 

Next difficult period has mainly focused on 
phenomenological manifestations of instanton- 
induced effects. The so called 'instanton liquid' 
model [|| has emerged as a qualitative picture. It 
suggested relative diluteness and large action per 
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instanton due to their relatively small size, but 
also emphasized significant interaction in the en- 
semble as the origin of density stabilization. It 
was shown, that this picture can explain several 
puzzles of the hadronic world (see e.g. @Jl[] for 
details). 

Attempts to describe interacting instantons 
were initiated by the variational approach ||], 
which has qualitatively reproduced the 'instan- 
ton liquid' picture. Further numerical studies 
of this problem have allowed to get rid of 
many approximations and eventually included 
fermionic effects to all orders in 't Hooft effective 
Lagrangian We return to discussion of this 
approach below, and now let me jump directly 
to several important steps made during the last 
year. 

First of all, more than 40 correlation functions 
were calculated in the framework of the simplest 
ensemble of the kind, the Random Instanton Liq- 
uid Model (RILM) [§. Agreement with data is 
generally good, and in some cases (including 7r,N 
etc) it is really surprising. 

Next, some of those functions were calculated 
on the lattice [|l5| , also with good agreement with 
the RILM results. That the agreement is not 
ocasional but based on the same dynamics was 
clearly shown by recent findings, a kind of deci- 
sive experiment, reported at this conference [l6| . 
'Cooled' lattice configurations, containing only 
instantons, not only has reproduced parameters 
of the 'instanton liquid', but they also lead to 
about the same correlation functions! 

I think these studies, taken together, has essen- 
tially answered many old questions like: why a 
nucleon is bound? Instanton-induced attraction 
rather than confinement or perturbative effects 
clearly play the major role here. 

2. WHY INSTANTONS? 

The main formal reason why instantons are so 
important for physics of light fermions is related 
to the famous 't Hooft 'zero modes', the localized 
solutions of the Dirac equation 

Df,'yf,(j)o{x) = 



in the instanton field. Evaluating the (Euclidean) 
quark propagator S = — l/[i£'^7^ -I- im] for m — > 
one has to deal mainly with small eigenvalues. 

It is convenient to look at the instanton as a 
'trap' for quarks, something like a 'receptor' atom 
in a semiconductor, creating a new state with the 
energy value being forbidden otherwize. At finite 
density of such atoms, an electron can propagate 
far, just by hopping from one atom to another. 
The same is true for finite instanton density, lead- 
ing to a 'zero mode zone' of collectivized quark 
states. That is the mechanism leading to the non- 
zero quark condensate. If more than one quark 
is travelling in the QCD vacuum {qq for mesons 
and qqq for baryons), they 'hop' over the same 
instantons. It implies an attractive interaction, 
which is in fact the one binding quarks together. 

Why instantons and not any other fluctuation 
of the gauge field? Here one has to look more 
specifically into the chiral and flavor structure of 
the instanton-induced interaction. As shown by 't 
Hooft, at tunneling quarks with one chirality 'dive 
into the Dirac sea' while those with the opposite 
chirality 'emerge' from it. Therefore instanton- 
induced forces should be stronger in scalar and 
pseudoscalar channels compared to vector or axial 
ones. Looking at phenomenological correlators at 
small distances, one finds that it is exactly right. 

Consider corrections to correlation functions 
are provided by instantons to the free propaga- 
tion at small distances. (If those corrections are 
relatively small, and this is the main point, one 
may hope to use the 't Hooft interaction in the 
lowest order, which makes consideration simple.) 
There are 4 such channels for 2 flavors and it is 
a simple matter to see that correction is positive 
(or attractive) for tt, a channels and negative (or 
repulsive for 6, rf ones. Thus the same mechanism 
leads to both light pion and heavy rj'l Both split- 
ting from 'typical mesons' are large, which is a 
very strong hint. 

3. CORRELATORS IN THE INSTAN- 
TON VACUUM 

Now we proceed from qualitative hints to quan- 
titative calculations. We have to evaluate a quark 
propagator in the multi-instanton field configura- 
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tion, which can be done as foUows: 

Of ^ V Mx)(t>t{y) , f ^^ 

b[x,y) = i^zMZ — ^ : \- iONZM(x,y)) 

A — im 

where the first term is the sum over states be- 
longing to the 'zero mode zone'. The non-zero 
modes (analogs of 'scattering states') are taken 
into account by the last term, see details in |^ . 

We have first calculated correlators for the sim- 
plest ensemble possible, the 'random instanton 
liquid model' (RILM), in which: (i) all instantons 
have the same size po = .35fm; (ii) they have 
random positions and orientations; (iii) instan- 
ton and anti-instanton densities are equal, and in 
sum it is no = 1/to~'*. These are the param- 
eters suggested a decade ago in Q, the density 
comes from the gluon condensate and size from 
various other things, say from the quark conden- 
sate value. 'Cooling' of lattice configurations pro- 
vide a way to see instantons and check them, see 
|l6| for the latest results and earlier references. 
Let me only state that these numbers are essen- 
tially confirmed. 

The main step in the calculation is inversion of 
the Dirac operator, written in the zero-mode sub- 
space. (We typically use in sum 256 instantons 
and anti-instantons, which tells the dimension of 
this matrix and the volume of the box.) 

Although the quark propagators are gauge de- 
pendent, we have looked at them first in order 
to see whether they can be reproduced by any 
simple model, say by 'constituent quarks' with a 
constant mass. We have found that chirality-non- 
flipping part of the propagator indeed looks as if 
quark get a mass about 300-400 MeV, but the 
chirality-flipping part does not look like that at 
all. None of many correlators calculated behave 
like 'constituent quark model'. 

Our results for tt, p, N, A channels are shown 
in Fig.|l| and Fig.|| All correlators are plotted in a 
normalized way, divided by those correcponding 
to free quark propagation: that is why all of them 
converge to 1 at small distances. Solid lines cor- 
respond to experiment , while the long-dashed 
and short-dashed curves correspond to QCD sum 
rule predictions and [|^, respectively. 

The agreement for the pion curve is as perfect 
as it can be: both the mass (142±12 MeV) and 



Figure 1. Correlators for tt, p channels according 
to RILM results (the closed triangles) and lattice 
calculation (open squares) versus distance x in 
fm. 



the (pseudoscalar) coupling are reproduced cor- 
rectly, inside the error bars! Note large devia- 
tions from perturbative behavior happens at very 
small distances for the tt channel, while exactly 
the opposite is observed in the p case, the plotted 
ratio remains close to 1 up to a very large x. Both 
RILM and lattice has reproduced that non-trivial 
observation. 

Proceeding to baryonic channels, let me men- 
tion that we have actually measured all 6 nucleon 
correlators and 4 delta ones, and have fitted them 
all. Again, agreement between RILM and lat- 
tice results is surprisingly good, literally inside 
the error bars. Both display a qualitative differ- 
ence between the nucleon and the delta correla- 
tors: this can be traced to attractive instanton- 
induces forces for the spin-isospin-zero diquarks 
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Figure 2. Same as in Fig.l but for N,A channels 



0. Without any one-gluon exchange the RILM 
predicts the — A sphtting (actually, we have 
found that in RILM rriN = 960 ± SOMeF and 
mA — 1440 ± 70MeV, so the splitting is in fact 
somewhat too large). 

There is no place here to discuss many other 
channels. The most difficult case proved to be 
the isosinglet scalar a, for which one should not 
only evaluate the 'double quark loop' term, but 
also subtract the disconnected | < qq > |^ part. 
Curiously enough, we have found dominance of a 
light state, with m~ 500 MeV, reminiscent of the 
'sigma meson' of 60's. For several reasons such 
measurements are now beyond the reach of lat- 



tice calculations, although existence of attractive 
interaction at a: 1/2 fm can probably be seen. 

Not included in the preprints mentiond are 
studies of the the so called 'wave functions' 
(known also as Bethe-Salpeter amlitudes). They 
have shown with greater clarity that RILM 
does lead to quark binding, without perturbative 
(Coulomb-like) and confining forces. Although 
the shape of the wave function is not the same as 
was found on the lattice, its width is only slightly 
larger. The main qualitative features (e.g.: n,N 
are more compact than p, A) are also reproduced. 

4. INTERACTING INSTANTONS 

Clearly, RILM discussed above cannot be but 
a crude approximation: at least, the very phe- 
nomenon studied above, 'hopping' of quarks 
from one instanton to another, should produce 
strong correlation between them. Another obvi- 
ous source of interaction is the non-linear gluonic 
Lagrangian: a superposition of instantons and 
anti-instantons have the action different from the 
sum of the actions. Effective statistical system 
should be described by a partition function 
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(where dfl is the measure in space of collective 
coordinates, 12 per instanton). It is a problem 
similar to those traditionally studied in statisti- 
cal mechanics, the compication is the fermionic 
determinant which is non-local. However, it is 
still orders and orders of magnitude simpler than 
the lattice gauge theory. As shown in 12| this 
statistical sum describes a liquid, in which chi- 
ral symmetry is broken. Recent studies of cor- 
relation functions with interacting ensemble have 
shown other significant improvement over RILM. 
In particular, the global fluctuations of the topo- 
logical charge are screened, and for m — > the 
topological susceptibility vanishes, as it should. 
The distribution of Dirac eigenvalues A has a dif- 
ferent shape, with a dip at A = instaed of a 
peak, and satisfy several non-trivial general the- 
orems, see pof . 

Unfortunately, it is still not quite quantita- 
tive theory. The reason is we do not know how 
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to separate 'semiclassicar fluctuations with large 
action, from very close instanton-anti-instanton 
pair on the bottom of the 'valley' ([^). Those 
configurations has small field, therefore they 
should be included in the perturbation theory. If 
we suppress them by some repulsive core, we get 
reasonable results, but a satisfactory solution of 
this problem is still missing. 

5. CHIRAL SYMMETRY RESTORA- 
TION 

Instantons at finite temperatures is a subject 
for separate talk, and here let me only discuss 
recent development related with critical phenom- 
ena near the chiral symmetry restoration point. 

The main phenomenon in this region is a strong 
'pairing' of instantons, leading to splitting of the 
'instanton liquid' into the // molecules. The first 
(strongly simplified) discussion of chiral restora- 
tion transition at this angle was made in fig] ]. 

One new finding is strong and rapid 'polar- 
ization' of these molecules in the critical re- 
gion. // interaction at finite T were stud- 
ied in details in |p7| , and the main anisotropy 
comes from the quark-induced interaction ~ 
\sin{TrTT)/cosh{TTTr)\^'^f , (where r,r are dis- 
tance between the two centers in time and space 
directions). Around T — c the point r = 0,r = 
1/(2T) is strongly prefered, and, as shown in 
Fig.^, the degree of polarization rapidly frows. 

This polarization leads to rapid growth of the 
gluonic energy density in the transition region, 
which was in fact observed in lattice simulations 
with dynamical fermions. Recall that, in terms 
of Minkowski field strengths, it is 
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, (ll/3)jV, - (2/3)jVj , , 
^ 1287r2 ^ 



For 'unpolarized' objects the first term is zero, 
and only the second ('anomalous') one con- 
tributes. However, for the 'polarized' molecules 
one finds E'^ = —.8B^, and the first term works. 
It produces a jump in e(T) of the order of 1 
GeV/fm^. 

The second interesting findings are specific 
changes in the correlation functions below Tc- A 



Figure 3. Fraction of the // molecules for which 
r > r for QCD (squares) and fir 4 massless quark 
flavors (stars). 



simple 'cocktail model' Q was used, contain- 
ing both 'random' component and some 'molecule 
fraction' with weight /,„ = 2Nmoiecuies/Nan. We 
have first found that < qq >, depends on fm in 
a way similar to its T-dependence, measured on 
the lattice: it changes little first, and then rapidly 
vanishes at /„ — > 1. Different correlation func- 
tions depend on /,„ quite differently. Two exam- 
ples are shown in Fig.^ for tt, p channels. One 
can see, that while the vector correlator is not 
changed much, the pion on display remarkable 
stability for /,„ = — 0.8, with subsequent strong 
drop toward /,„ = 1. At the last point com- 
plete chiral symmetry get restored, so the pion 
correlator coinsides with its scalar partners cr, 6. 
More studies along these lines and their compar- 
ison with recent lattice works are in progress. 
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